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Abstract 

We address a simple connection between results of Hamiltonian nonlin- 
ear dynamical theory and thermostatistics. Using a properly denned dy- 
namical temperature in low-dimensional symplectic maps, we display and 
characterize long-standing quasi-stationary states that eventually cross 
over to a Boltzmann-Gibbs-like regime. As time evolves, the geometrical 
properties (e.g., fractal dimension) of the phase space change sensibly, and 
the duration of the anomalous regime diverges with decreasing chaoticity. 
The scenario that emerges is consistent with the nonextensive statistical 
mechanics one. 

PACS numbers: 05.70.Ln 05.10.-a 05.20.Gg 05.45.Ac 

The methods of usual, Boltzmann-Gibbs (BG), statistical mechanics apply 
to impressively large classes of macroscopic systems. However, the situation 
is more delicate for complex systems. Indeed, turbulent fluids |XJ, high-energy 
collision processes j2i , classical |Sj and quantum chaos @J , stellar self-gravitating 
systems jS], granular systems UJ, economics 0, motion of micro-organisms 
jS], and others, frequently exhibit anomalous behaviors where alternative ap- 
proaches arc needed. In particular, in many-body long-range- interacting Hamil- 
tonian systems, it has been recently observed the emergence of long-standing 
(in the thermodynamical limit infinite-lasting) quasi- stationary (metastable) 
states (QSS) characterized by non-Gaussian velocity distributions, before the 
Boltzmann-Gibbs (BG) equilibrium is attained [HI EI]. This is a major concern, 
as, for these Hamiltonian systems, the foundation of the BG equilibrium ther- 
modynamics is questioned. Using standard results, in this letter we address a 
simple connection between chaos theory and thermostatistics, and we focus on 
a paradigmatic dynamical mechanism that produces QSS very similar to those 
detected in EJ. These QSS are displayed and characterized by means of 
low-dimensional symplectic maps. 

*E-mail addresses: baldovin@cbpf.br, edgardo@cbpf.br, tsallis@cbpf.br 



The foundation of the Boltzmann-Gibbs (BG) equilibrium thermodynam- 
ics lies on a sufficiently complete and uniform occupation of the system phase 
space (the finite Lebesgue measure T-space), taking into account symmetry, en- 
ergy and similar restrictions. The BG equilibrium descends in fact from the 
equal- a-priori-probability postulate, that characterizes the microcanonical Gibb- 
sian ensemble (see, e.g., According to this postulate, each equally-sized 

accessible region of the phase space (under the macroscopic conditions of the 
system) equally likely contains the microscopic state of the system. As Einstein 
pointed out in his criticism of the Boltzmann principle S = klnW this 
postulate should not be taken a priori, but rather justified a posteriori by the 
underlying dynamics. Indeed, if dynamics is sufficiently chaotic, large portion 
of the phase space are rapidly occupied by the trajectory of the system and the 
postulate is a very accurate representation of the dynamical behavior, as testi- 
fies more than a century of successes of the BG formalism. But there are also 
many situations where the system displays an intricate dynamical behavior, as 
it happens for example at the border between regular and chaotic regimes. 

At this border, for a large class of Hamiltonian systems, a mechanism based 
on the KAM theory operates, which we briefly review now. A continuous Hamil- 
tonian system with n degrees of freedom may be written in the form 

H = H (h, ...,/„) + e V(h,6 u I n ,0 n ), (1) 

where H is integrable (ii, ...,/„ are its integrals of motion), e << 1, and V is 
a nonlinear perturbation. Under certain hypothesis (see, e.g., |T3]). for e = 
the trajectories lie on invariant n-dimensional tori. A special subset of these 
tori are called resonance tori. Specifically, if we introduce the (non-degenerate) 
frequencies of the unperturbed motion: uij = ^j - (j = 1, ...,n), we have that 

the condition Y^j=i m j UJ j ~ (where m,j are integer numbers) defines the res- 
onance tori. Each resonance torus involves the formation of a separatrix loop. 
The action of the perturbation, for small enough e ^ 0, deforms normal tori 
into KAM-tori, and, in correspondence with the resonance tori, destroys the 
separatrices replacing them with stochastic layers. Resonance tori, in the space 
spanned by u>±, ..,w n , lie in the intersection between the hyperplane defined by 
the resonance condition and the hypersurface of energy E = Hq{ui\, ...,u> n ). In 
the case n > 2, resonance tori must, for topological reasons, intersect between 
them. Consequently, while for n < 2 the stochastic layers are distinct for e suf- 
ficiently small, for n > 2 they merge into a single connected stochastic web that 
is dense in the phase space for all e / 0, and there is room for Arnold diffusion 
processes. We remark that, for n = 2, KAM-tori constitute total barriers for 
diffusive processes in the phase space; nevertheless, inside the stochastic sea, it 
is possible to find Cantor sets, named cantori, that constitute partial barriers 
for diffusion (see ^3] for details). 

A convenient way of studying Hamiltonian systems is by using symplectic 
maps. A (2n — 2)-dimensional symplectic map is obtained from conservative 
Hamiltonian systems with n degrees of freedom by taking a Poincare section 
over the hypersurface of constant energy. Interestingly enough, a 2n-dimensional 
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symplectic map is also the result of a Poincare section on the phase space of an 
open system of n degrees of freedom with a Hamiltonian that depends periodi- 
cally on time. We remark that in both cases the map has a symplectic structure; 
this assures (hyper)volume conservation in the phase space. The advantage of 
maps lies on the reduced dimension of the phase space and on the use of a 
discrete time. In this letter we specifically address some symplectic (hence con- 
servative) maps in order to discuss how equilibrium and quasi-equilibrium can 
be attained in phase space. 

Let us start through the analysis of a prototypical 2-dimcnsional symplectic 
map (n < 2), the standard (or kicked rotor) map 



(a G R, t — 0, 1, ...). 2irp may be regarded as an angular momentum variable. 
Notice that, consistently with our scope, we have used the symmetry properties 
of the map and defined, as usually, the angular momentum mod 1. The standard 
map is integrablc for a = 0, while chaoticity rapidly increases with \a\. 

In [H1E3|j the emergence of the dynamical QSS appeared to be dependent 
on the initial conditions. Specifically, it was shown that for some classical long- 
range-interacting iV-rotor Hamiltonian models, a basin of attraction of initial 
data exist for which the system dynamically evolves into a QSS whose duration 
diverges as N — > oo. Typical examples of this basin of attraction are out-of- 
equilibrium initial conditions called 'water bag' initial conditions, characterized 
by a uniform initial distribution of the angular momenta around zero (see [9"l ll()| 
for details). In the case of the standard map, we first observe that the points 
(0, 1/2) and (1/2, 1/2) are a 2-cycle for all a so that we can use them as refer- 
ential for studying the properties of the phase space with respect to variation of 
the parameter a. With some analogy with [91110). our out-of-equilibrium 'water 
bag' initial conditions are defined by considering at t = a statistical ensemble 
of M copies of the standard map with arbitrary 6 and p randomly distributed in 
a small region around p = 1/2. In standard statistical mechanics, when dealing 
with systems with diagonal kinetic matrix and zero average momentum, the 
temperature is proportional to the average square momentum per particle. As 
we analyze situations with nonzero 'bulk' motion, the analogous concept, which 
we shall refer to as (dimensionless) 'dynamical temperature', can be defined as 
the variance of the angular momentum: T = ((p— (p)) 2 ) = (p 2 ) ~ (p) 2 : where () 
means ensemble average. The temperature associated with the uniform ensem- 
ble (that we will call BG temperature because of its similarity with the equal- 

a-priori-probability postulate) is given by Tbg = Jq dp p 2 — ( dp p^j = 1/12. 

It should be noticed that in the present conservative model, the 'temperature' 
T is necessarily bounded since p itself is bounded, in contrast with a true ther- 
modynamical temperature, which is of course unbounded. For large values of 
\a\ (i.e., strong chaoticity) in map J5J), the temperature of the 'water bag' initial 
ensemble rapidly relaxes to Tbg (let us stress that the subindex BG stands for 



9{t + l) 



p(t)-i sin[27T0(t)] + 6(t) (modi), 



(2) 



p(* + l) 
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the fact that it corresponds to uniform occupation of the accessible phase space 
of the map, to be distinguished from the phase space of the physical kicked 
rotor from which the map was originally deduced). Our aim is to study what 
happens in the transition to regularity obtained reducing the value of \a\ to- 
wards a = 0. In Fig. ^a) we see that the first effect of the reduction of \a\ is 
that limj^oo T(t) < Tbg- This is easily understood as follows. As chaoticity 
reduces, total barriers (KAM-tori) appear in the phase space. As a consequence, 
the points of the ensemble are prevented to reach all the regions of the phase 
space and the projection of the ensemble on the p axis produces a probability 
distribution function (PDF) with a variance smaller than the one of the uniform 
distribution. For values of a of order a ~ 1, a QSS emerges before the relax- 
ation to the final temperature. In fact, inside the stochastic sea partial barriers 
(cantori) begins to appear, and the initial 'water bag' first rapidly diffuse inside 
an area delimited by cantori and then slowly crosses over to the final relax- 
ation temperature. Fig. Qfc) illustrates this behavior. To obtain a quantitative 
description of the relaxation time we have plotted the iteration time in a loga- 
rithmic scale. We define the crossover time t c as the inflection point of the curve 
and we observe that it diverges, as a tends to a c = 0.971635406... from above, 
like t c ~ l/(a — a c ) 2 ' 7 (see Fig. ^b)). Just below this critical value in fact, 
the strongest cantori close ^2], and the relaxation to a higher temperature is 
prevented. Reducing a to smaller values causes the formation of more and more 
total barriers, so that the 'dynamical temperature' tends to zero for \a\ — > 0. 
^From the mechanism we have displayed, it is clear that it is possible to obtain 
these types of QSS even with other sets of initial conditions. Typically, it is 
sufficient to have the initial data localized inside the first partial barriers. In 
other words there is an entire basin of attraction of out-of-equilibrium initial 
conditions that leads to the formation of a certain kind of QSS. 

As we pointed out previously, the topology of the phase space changes dra- 
matically for n > 2. To address this case, we move next to a 4-dimensional 
symplectic map composed by two coupled standard maps: 



where a\,a 2 ,b 6 R, t = 0,1,..., and all variables are defined mod 1. If the 
coupling constant b vanishes the two standard maps decouple; if b = 2 the 
points (0, 1/2,0, 1/2) and (1/2,1/2, 1/2,1/2) are a 2-cycle for all (ai,o 2 ), hence 
we preserve in phase space the same referential that we had for a single standard 
map. For a generic value of b, all relevant present results remain qualitatively 
the same. Also, we set a\ = a 2 = a so that the system is invariant under 
permutation 1 <-> 2. Since we have two rotors now, the 'dynamical temperature' 
is naturally given by T = \ (< p\ > + < p\ > ~ < pi > 2 ~ < p 2 > 2 ), hence the 
BG temperature remains Tbg = 1/12. 



Pi(* + 1) 
2 (* + l) 
P 2 (i + 1) 



Pl (t + l) + 9 1 (t) + bp 2 (t + l), 
Pi(t) + ^wn[2n9i(t)], 

p 2 (t + i) + e 2 (t) + b Pl (t + i), 

p 2 (t) + ^sm[2ne 2 (t)l 



(3) 
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As before we consider 'water bag' initial conditions, i.e., an ensemble of M 
points with arbitrary (9\, O2), and angular momenta randomly distributed inside 
a small region around p\ = P2 — 1/2. The result is qualitatively similar to the 
one displayed in Fig. ^a) for a > a c . Large values of a correspond to Tbg and 
reducing a we observe the formation of a QSS that, after some time, relaxes to 
a temperature T < Tbg- The first major difference with the 2-dimensional case 
is that, because of the Arnold diffusion processes, the relaxation to a higher 
temperature occurs (waiting enough time) for all a ^ (i.e., a c ~ 0, defining a c 
as the value where t c diverges). Moreover, the reason why the ensemble does not 
relax to the BG temperature for small values of a is here quite different from that 
for the 2-dimcnsional case. Indeed, with this choice of initial data, the initial 
'water bag' intersects at least a macroscopic island, as clearly appreciated in Fig. 
IH(al); the points that are set inside the island do not diffuse to the outside. As 
a result, the projection of the ensemble on the plane p\,P2 conserve a denser 
central part for all times (Fig. El(bl), (cl)). 

If we instead shift the initial 'water bag', say towards the lower part of the 
phase space, we can set the points outside this island (Fig. |21 (a2)-(c2)). In this 
case we obtain a crucial qualitatively new phenomenon, namely the formation, 
for small values of a, of a QSS that eventually relaxes to the BG temperature 
(see Fig. Eta)). Notice also that the crossover time t c diverges, as a — > 0, faster 
than in the 2-dimensional case: t c ~ 1/a 5,2 (see Fig. ED 5 ))- We remark that 
the relaxation to the BG temperature occurs here even if the presence of islands 
in the phase space violate the equal a priori postulate. In other words, it is 
possible to obtain a weak violation of the postulate that does preserve a uniform 
distribution once the ensemble is projected over the plane pi,P2, in the same 
sense that a sponge projects a uniform shadow on a wall. These QSS can in 
fact be geometrically characterized by the fractal dimension df. Overcoming 
some numerical difficulties involved in a fractal analysis in 4 dimensions, we 
illustrate what happens in Fig. Etc), constructed using a box-counting algorithm 
[T5] (in fact, the phase space exhibits strong inhomogeneities which suggest a 
multifractal structure). During the QSS the ensemble is, for small a, associated 
with a nontrivial fractal dimension df ~ 2.7, while, once it crosses over to 
the BG-like regime, it distributes itself occupying the full dimensionality of the 
phase space, thus attaining df = 4. 

These simple models allows also for the discussion of different types of QSS. 
For instance, at t = we can set up 'double water bag' initial conditions consid- 
ering an ensemble of M copies of two coupled standard maps with arbitrary 
61,62 and angular momenta randomly distributed inside two small regions: 
Pi,P2 = + S and pi,p2 = 1 — 8 (0 < S « 1). In this case, the initial 
temperature T(0) is higher than Tbg, because the PDFs projected on the Pi,P2 
axes are double peaked. Relaxation to Tbg occurs then from above, as can be 
seen in Fig. This is precisely in what the phenomenon observed in [2J differs 
from the one in UJ. 

It is relevant to notice that, as a consequence of Arnold diffusion, our results 
for the 4-dimensional map capture the qualitative behavior of higher-dimension 
symplcctic maps. In fact, the two-plateaux structure has been confirmed by 
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numerical analyses where hundreds of standard maps arc all-with-all coupled as 
in Eq. © [SI. 

Through the numerical analysis of low-dimensional symplectic maps, we 
have displayed how complex paradigmatic structures associated with conser- 
vative nonlinear dynamics can generate anomalous thermodynamical behavior. 
Particularly, we have exhibited and studied the emergence, while approaching 
integrability (i.e., when chaoticity decreases), of QSS suggestively similar to 
those observed in long-range iV-body systems (a playing a role analo- 
gous to 1/N). A central result is that, in contrast with what happens for the 
BG equilibrium, these QSS correspond to a nontrivial fractal dimension. This 
situation reminds the phase space structure of logistic-like maps at the edge of 
chaos (also characterized by a nontrivial fractal dimension), where exact ana- 
lytical connections with the nonextensive statistical mechanics |17| have been 
established [jjj. 
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Figure 1: (a) Time evolution of the dynamical temperature T of a standard map, 
for typical values of a. We start with 'water bag' initial conditions (M — 2500 points 
in < 8 < 1, p = 0.5 ± 5 10~ 4 ). In order to eliminate cyclical fluctuations, the 
dots represent average of 10 iteration steps; moreover, each curve is the average of 50 
realizations, (b) Inverse crossover time t c (inflection point between the QSS and the 
BG regimes) vs. l/(a — a c ) 2,7 . No inflection points subsist if t is linearly represented, 
(c) Time evolution of the the ensemble in (a) for a = 1.1 (first row) and PDF of its 
angular momentum (second row), t = 0: 'water bag' initial conditions; t = t\ = 500: 
the ensemble is mostly restricted by cantori; t = ti = 10 5 : the ensemble is confined 
inside KAM-tori. 
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Figure 2: Phase space analysis of the evolution of 'water bag' ensembles for two 
coupled standard maps for (a,b) = (0.4,2). First row: 'Water bag' initial conditions 
< 61,82 < 1, Pi,P2 = 0.5 ± 5 10~ 3 . Second row: 'Water bag' initial conditions 
< 61,82 < 1, Pi,P2 = 0.25 ± 5 10~ 3 . (a) Projection on the (81, J5i)-plane of the 
central slice of the phase space (02, P2 = 0.5 ± 10~ 2 ), for the orbit < t < ti = 10 4 . 
(c),(c) Projection on the (pi,p2)-plane ofgwhole phase space for the iterate at time 
t 2 = 15 and t 3 = 2 10 4 . 
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Figure 3: (a) Time evolution of the dynamical temperature T of two coupled standard 
maps, for b — 2 and typical values of 5. We start with water bag initial conditions 
(M = 1296 points with < 1 ,0 2 < 1, and p\,pi = 0.25 ± 5 10" 3 ); moreover, an 
average was taken over 35 realizations. See Fig. [5] for ti and tz. (b) Inverse crossover 
time t c vs. 1/5 5,2 . (c) Time evolution of the fractal dimension of a single initial 
ensemble in the same setup of (a). 
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Figure 4: Same as Fig. l3Ta).(b) but with 'double water bag' initial conditions: 
< 01,02 < 1; Pi,P2 randomly distributed inside one of the two regions pi,p2 = 
+ 10~ 2 ,pi,p 2 = l-10- 2 . 
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